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1 Introduction 

This work is concerned with existence and uniqueness of solutions to stochas- 
tic porous media equations 

dX(t) - AV(X(t))dt = B(X(t))dW(t) in (0, T) x ■= Q T , 
V(X(t)) =0 on (0, T) x dG := S T , (1.1) 
X(0) = x inff, 

where & is an open, bounded domain of R d , d > 1, with smooth boundary 
dG, W(t) is a cylindrical Wiener process on L 2 (G), B : H -> L(L 2 (G), L 2 (<ff)) 
is a Lipschitz continuous operator to be precised below and H : = 
The function \P : R — > R (or more generally the multivalued function 
: R — > 2 M ) is a maximal monotone graph in R x R. 

Existence results for equation (11. ip were obtained in [8] (see also [3],[1]) 
in the special case B = \fQ, with Q linear nonnegative, Tr Q < +oo and 
\1/ G C 1 (R) satifying the growth condition 

h + hlsl"- 1 < V'(s) < A; 2 (l + |s| r — 1 ), sGR, (1.2) 

where ki, k 2 > 0, /c 3 G R, r > 1. 

Under these growth conditions on equation (11.11) covers many impor- 
tant models of dynamics of an ideal gas in a porous medium (see e.g. [T]) 
but excludes, however, other significant physical models such as plasma fast 
diffusion ([5]) which arises for \I/ (s) = y/s, phase transitions or dynamics of 
saturated underground water flows (Richard's equation). In the later case 
multivalued monotone graphs \I/ might appear (see [I2])- Recently in [15] 
(see also [2]) the existence results of [8] were extended to the case of mono- 
tone nonlinearities \l/ such that s h- > ^(s) is (comparable to) a A 2 -regular 
Young function (cf. assumption (Al) in [15]) thus including the fast diffu- 
sion model. As a matter of fact, in the line of the classical work of N. Krylov 
and B. Rozovskii [10] the approach used in [15] is a variational one i.e. one 
considers the stochastic equation (11.11) in a duality setting induced by a func- 
tional triplet V C H C V and this requires to find appropriate spaces V and 
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H . This was done in [T5] in an elaborate way even with A replaced by very 
general (not necessarily differential) operators L. 

The method we use here is quite different and essentially an L 1 -approach 
relying on weak compacteness techniques in L 1 (Qt) via the Dunford- Pettis 
theorem which involve minimal growth assumptions on ^ . Restricted to sin- 
gle valued continuous functions \& the main result, Theorem 12.21 below, gives 
existence and uniqueness of solutions only assuming that hm s _^ +0O ^f(s) = 
+00, lims^-oo = —00, \1> monotonically increasing and 

hmsup fs lW ,N < +°°- I 1 - 3 ) 
M-H-00 Jo V{t)dt 



We note that the assumptions on \1/ in [15]) imply our assumptions. In this 
sense, under assumption (H2) below in the noise, the results on this paper 
extend those in [15] in case L = A if G is bounded and if the coefficients do 
not depend on (t,u). The latter two were not assumed in [15]. On the other 
hand a growth condition on \l/ is imposed in [15] (cf. [T5l Lemma 3.2]) which 
is not done here. Another main progress of this paper is that \1/ is no longer 
assumed to be continuous, it might be multivalued and with exponential 
growth to ±00 (for instance of the form exp (a|x| p )). We note that (11.31) is 
not a growth condition at +00 but a kind of symmetry condition about the 
behaviour of \& at ±00. If \& is a maximal monotone graph with potential j 
(i.e. = dj) then ( 11.3j) takes the form (see Hypothesis {H3) below) 

hmsup < +00. 

|s|->+oo ]{ S ) 

Anyway this condition is automatically satified for even monotonically in- 
creasing functions \P or e.g. if a condition of the form (11.21) is satisfied. We 
note, however, that because of our very general conditions on \1/ the solution 
of (11. ip will be pathwise only weakly continuous in H. It seems impossible 
to prove strong continuity. 

1.1 Notations 

& is a bounded open subset of R d , d > 1 with smooth boundary d& . We set 
Q T = (0, T) x 0, S T = (0, T) x dff. 
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Moreover L p (0), L p (Qt), p > 1, are standard LP- function spaces and Hq(&), 
H k (@) are Sobolev spaces on 0. Denote by H := H~ x {0) the dual of H\{G) 
with the norm and the scalar product 

:= {A~ l u, u) 1/2 , (u, = {A~ x u, v), 

where (-, ■) is the pairing between Hq(0) and Hq 1 (0) which coincides with 
the scalar product of L 2 (<ff). We have denoted by A the Laplace operator 
with Dirichlet homogeneous boundary conditions, i.e. 

Au = -Au, ueD(A) = H 2 (0)nH^(ff). (1.4) 

Given a Hilbert space U, the norm of U will be denoted by \ ■ \u and the 
scalar product by (-,•)*/■ By C([0,T];U) we shall denote the space of U- 
valued continuous functions on [0, T] and by C w ([0, T]; U) the space of weakly 
continuous functions from [0, T] to U . 

Given two Hilbert spaces U and V we shall denote by L(U, V) the space 
of linear continuous operators from U to V and by Lhs(U, V) the space of 
Hilbert-Schmidt operators F : U — > V with the norm 

(oo \ V 2 

where {e;} is an orthonormal basis in U. 

If j : R — > (— oo, +oo] is a lower semicontinuous convex function we 
denote by <9j : R — ► 2 K the sub differential of j, i.e. 

dj(y) = {9eR: j(y)< j{z) + 9{y - z), V z el} 

and by j* the conjugate of j (the Legendre transform of j), 

j*(p) = snp{py-j(y) : y £ R}. 

We recall that <9j* = (dj)" 1 (see e.g. [2], [6]), 

j{y) + j*(p) = PJ/ if and only if p e dj(y) (1.6) 

and 

j(it) + j*Qo) > for all p, u G R. (1.7) 
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Moreover :— dj is maximal monotone, i.e. 

(yi - V2)(pi - P2) > for all pi G dj(yi), i = 1, 2 

and R(l + dj) = R. 

Given a multivalued function $ : R — > 2 K we shall denote by -D($) = {u£ 
R : 7^ 0} the domain of $ and by = {f : v E u G -D($)} 

its range. 

Given a maximal monotone graph \P : R — > 2 R there is a unique lower 
semi continuous convex function j : R — > (— 00, +00] such that ^ := <9j. The 
function j is unique up to an additive constant and called the potential of 

For any maximal monotone graph \1/ and A > by 

* A = 7 (1 - (1 + A^r 1 ) G $(1 + A*)" 1 
A 

we denote the Yosida approximation of Here 1 stands for the identity 
function. ty\ is Lipschitzian and monotonically increasing. 

2 The main result 

2.1 Hypotheses 

(Hi) W(t)is a cylindrical Wiener process on L 2 (ff) defined by 

00 

W(t) = J2Pk(t)e k , (2.1) 

k=l 

where {fik} is a sequence of mutually independent Brownian motions on a 
filtered probability space (Q, J^", {^t}t>o, 1P)> with right continuous filtration 
and {ek} is an orthonormal basis in L 2 (ff). To be more specific {e k } will be 
chosen as the normalized sequence of eigenf unctions of the operator A, hence 
e k G L v (0) for allkeN,p> 1. 

(H 2 ) B is Lipschitzian from H = H~ x (&) to L HS (L 2 (ff), D(A^)) where 7 > 
d/2. 

(H3) ^ : R — > 2 R is a maximal monotone graph on R x R such that G \l/(0), 

D(V)=R, R(V) = R (2.2) 
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and 

i(-s) 

limsup ^— - < +00. (2.3) 

M-+QO j{ S ) 

Here j : R — > R is the potential of \P i.e. <9j = which under assumption 
( 12.21) is a continuous convex function. Since G ^(0), by definition we have 
j(0) = inf j. Hence subtracting j(0) we can take j such that j(0) = and 
j > 0, hence j* > j*(0) = 0. It should be recalled (see e.g. [2], [6]) that the 
condition R(^) = R is equivalent to 

j(y) < 00 V y G R, lim = +00 (2.4) 

bHoo 1 2/ j 

and that the condition D(^/) = R is equivalent to 

j*(y) < 00 V y G R, lim = +00 (2.5) 

I1/H00 \y\ 

Hypothesis (-£^3) automatically holds if ^ is a monotonically increasing, con- 
tinuous function on R satisfying condition (11.31) and 

lim ^{s) = +00, lim \I/ (s) = —00. 

In particular, it is satisfied by functions \l/ satisfying (jl.2p for r > or more 
generally by those satisfying assumption (Al) in [T5] . 

We need some more notations. Given a Banach space Z we shall denote 

by 

CV([0, T]; Z) = C([0, T]; L 2 (fi, J?, P; Z)) 

the space of all continuous adapted stochastic processes which are mean 
square continuous. The space 

L 2 W ([0, T};Z) = L 2 ([0, T];L 2 (Q, P : Z)) 

is similarly defined (see e.g. [7J, [H]). 

Definition 2.1 An adapted process X G C w ([0, T]; if) nL 1 ((0, T) x ^xfi) ; 
swc/i t/iat X G C w ([0, T], H),¥-a.s., is said to be a strong solution to equation 
(11.11) if there exists a process 77 G L 1 ((0,T) x G x f2) stzc/i £/iai 

7j(t,0 G a.e. (t,£) G Q T , P-a.s. (2.6) 
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['rj^dse C™{[Q,T]-Hl{0)), (2.7) 
Jo 

X{t)-A [ 7]{s)ds = x+ [ B(X(s))dW(s), V t G [0, T], F-a.s. (2.8) 
Jo Jo 

j(X),f( V ) G ^((O^x^xfi). (2.9) 

(Here J^rj(s)ds is initially defined as on L x (&)-valued Bochner integral). Of 
course, if ^ is single valued (12.61) - f)2.8p reduce to 

f m{X{s))ds G C w ([0, T]; if V)), (2.10) 

and 

X(t)-A [ V(X(s))ds = x+ f B(X(s))dW(s), V * G [0, T], P-a.s. 
Jo Jo 

(2.11) 

We note that X is as in Definition 12.11 is automatically predictable. 
Theorem 12.21 below is the main result of this work. 

Theorem 2.2 Under Hypotheses (Hi), (H2), {H3), for each x G H there 
is a unique strong solution X = X(t,x) to equation (11. ip . Moreover, the 
following estimate holds 

E\X(t,x)-X(t,y)\ 2 L 1 <C\x-y\ 2 _ 1 , forallt>0, (2.12) 

where C is independent of x,y G H . 

Theorem 12 . 21 will be proved in section 4 via fixed point arguments. Previously, 
we shall establish in section 3 the existence of solutions for the equation 

dY[t) - AV(Y(t))dt = G(t)dW(t) in Q T , 

(y(t)) = on E T , (2.13) 
Y(0) = x in G, 

where G : [0, T] — > Lhs(L 2 (&), D(A~')) is a predictable process such that 

E ||G(f)ll2 HS (L 2 (^),D(yl7))^ < +OO (2-14) 
JO 
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and 7 > d/2. By GdW we mean of course 

oo 

GdW = J2 Ge kdPk- 

k=i 

By a solution of (I2.13P we shall mean an adapted process Y satisfying along 
with 7] e L l ((Q,T) x x Q) conditions (j2T6>(j2^|) where B(X) is replaced 
byG. 

Theorem 2.3 Under Hypotheses (Hi), (H 3 ), (12.141) . for each x E H there 
is a unique strong solution Y = Ya(t,x) to equation (12 . 1 3[) in the sense of 
Definition \2.1\ Moreover, the following estimate holds 

E\Y Gl (t,x) - Y G2 {t,y)\\ <\x - y\\ 

r t (2.15) 

+E / \\G 1 (s)-G 2 (s)\\l Hs{LH ^ H) ds J for allt > 0, 
Jo 

for all x,y € H and G\, G 2 satisfying (12.141) . 

Remark 2.4 It should be noted that assumption (H2) excludes the case of 
equation (II. ip with covariance operator B of the form B(X) = X i.e. the 
case of multiplicative noise. However such an equation can be approximated 
taking B(X) = X * p € (p is a mollifier in G) or B(X) = (1 + eA)~ 5 X, e > 0. 

Remark 2.5 Assumption (^3) for example allows monotonically increasing 
functions ^ which are continuous from the right on M and have a finite 
number of jumps r 1; r 2 , ...,r^. However in this case one must fill the jumps 
by replacing the function \1/ by the maximal monotone (multivalued ) graph 
v]/(r) = ^(r) for r different from and \&(rj) = [W(rj) — ^(rj-i — 0)]. Such 
a situation might arise in modelling of underground water flows (see e.g. 
[T2"]). In this case ^ is the diffusivity function and (II. ip reduces to Richard's 
equation. It must be also said that Theorems 12.21 and 12.31 have natural 
extensions to equations of the form 

dX[t) - M!(X(t))dt + $(X(t))dt = B(X(t))dW(t), (2.16) 

where $ is a suitable monotonically increasing and continuous function (see 
|15j). Also as in [15] one might consider the case where \l/ = *&(X, u), u e fl, 
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but we do not go into details, here. We also note that assumption D(fy) = R 
in Hypothesis (H 3 ) excludes a situation of the following type 



Wfg) = I ^ for s < s °' = (°> +°°)' ( 2 17) 

U | = for s > s , ' 

where \& is a continuous monotonically increasing function such that oo) 
= — oo. In this case problem (11 .ip reduces to a stochastic variational inequal- 
ity and it is relevant in the description of saturation processes in infiltration. 
An analysis similar to that to be developped below shows that in this case in 
Definition 12.11 the solution is no more an L 1 - function but a bounded measure 
on Qt- We expect to give details in a later paper. 

Another situation of interest covered by our assumptions (see also [T5] ) 
is that of logarithmic diffusion equations arising in plasma physics see e.g. 
[13]. In this case \l/(s) = log(/x + \s\) sign(s). 



3 Proof of Theorem 2.3 



For every A > consider the approximating equation 

dX x (t) - A(^ x (X x (t)) + XX x (t))dt = G(t)dW(t) in (0, T) x := Q T , 
■V x (X x (t)) + \X x (t)) = on(0,T)x^, 
X x (0) = x in G, 

(3.1) 

which has a unique solution X x G CV([0,T]; H) such that 

X X ,V X (X X ) eL 2 w (O,T;H*(0)). 

Indeed, setting y x {t) = X x (t)-W G (t) where W G {t) = £ G(s)dW(s), we may 
rewrite (13. ip random equation 

y^(t) - AV x (y x (t) +W G (t)) = P-a.s. in Q T , 
*x(yx(t) + W G (t)) = on (0, T) x 80, (3.2) 
y x (0) =x in 0, 

where * A (y) = * A (y) + Xy, A > 0. Note that * A (0) = 0. 

For each u G the operator r(t) : Fq(^) -> H~ 1 {0) defined by 

r(t)j/ = -AW A ( 2 / + W G (t)), 
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is continuous, monotone and coercive, i.e. 

(T(t)y,y) > \\y+W G {t)\l l{e) -{V{t)y,W G {t)) > ± \y\l< {ff) -C x \W G (t)\ 2 H ^ y 

Then by classical existence theory for nonlinear equations (see e.g. 
equation (I3.2p has a unique solution 

y x G C([0, T];L 2 (0)) H L 2 (0, T; H 1 ^)) 

with ^ G L 2 (0, T; H~ l (ff)). The function X A (t) = y A (i) + W G (i) is of course 
an adapted process because the solution y\ to equation ( 13.21) is a continuous 
function of W G and so it satifies the requested condition. 

3.1 A-priori estimates 

From now on we shall fix uo G f2 and work with the corresponding solution 
yx to 1ET2]) . We have 

5 IvaWI-i + (*a(Va(*) + W o (*)) f y A (0 + Wfc(f)) 

^ dt (3.3) 

= (*x(yx(t) + W G (t)),W G (t)), 
which is equivalent to 
1 d 



2 dt 



\yx{t)\\ + (Vx(yx(t) + W G (t)),y x (t) + W G (t)) 

(3.4) 

-X(yx(t),yx(t) + W G {t)) + (*aM*) + Wg(*)), W G (t)). 



Now set Ja(w) = J^^x( r )dr and denote by j A the conjugate of j X - By (11.61) 
we have 

j^ x ( yx (t)+W G (t)))+jx(yx(t)+W G (t)) = *x(yx(t)+W G (t))(y x (t)+W G (t)). 
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Substituting this identity into (13.41) yields 

l\Vx(t)ti+ f [ (jx(yx(s) + W G (s))+f x (y x (y x (s) + W G (s)))dtds 
A Jo Jg 



- \x\ 
2 1 1 




, , (*x(y x {s) + W G (s))W G (s))d£ds 
'o Jg 



-A / / y x (s)(y x (s)+W G (s))d£ds, 
Jo J G 

Since j x is the Moreau approximation of j 



j x (u) = min jj'O) + — \u - v\ 2 : v G # j , 
we have (recall that the minimum is attained at v = (1 + A\1/) _1 m) 
j x (u) = + A*)- 1 ^) + \u - (1 + A*)"^! 2 , m 6 1. 



(3.5) 



We now set 



z x = (1 + A*)" 1 ^ + W G ), = *a(?/a + W G ). 



(3.6) 



(3.7) 



Then, using 1I3.6P and the fact that jl > j* for all A > 0, we see by ( 13. 5 p that 



^\yx{t)?-x 




(j(z x {s))+f(r] X ( S ))dZds 



JG 



+i k Jo L {yx{s) + Wg{s) ~ z ^ 2 ^ ds 



< 



+ [ [ V x(s)W G (s)d^ds - A / / y x (s)(y x (s)+W G (s))dtds. 
* Jo Jg Jo Jg 

(3.8) 

We now estimate the first integral from the right hand side of 1 13. 8ft as follows 




r} X (s)W G (s)d^ds 



JG 



<S [ [ \7} X (sMds, 

Jo Jg 



(3.9) 
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where 5 : = sup se r 0T i \Wg(s)\l°°(&) < +00. We note that by assumption 
(I2.14p and since 7 > d/2 it follows by Sobolev embedding that Wq(-) has 
continuous sample paths in D(A"<) C L°°(l?) and so 5 is indeed finite. 
Substituting in ($3M) yields 

J \yx(t)\\ + [ [ (j(z x (s))+f( Vx (s)dtds 
A Jo Jff 



1 



+ ^ / / (yx(s) + W G - z x (s)) 2 d^ds 



'0 </<^ 

"t 



< 

Since 



l\x\ 2 -! + S f f \vx(s)\d£ds-\ f [ y x (s)(y x (s) + W G (s))d^ds. 
1 Jo Je Jo Je 



-y x (s)(yx(s) + W G (s)) < -i |y A ( S )| 2 + ~ Wg( a ), 



we find 
1 
2 



f(j(z x (s))+f(r ]x (s))dtds+± f f \y x (s)\ 2 d^ds 
Jo Jff 1 Jo Je 



+ f / (!/A(a) + W G (s) - z x (s)) 2 d^ds 




Je 



2A 



1 a ■ 1" r a 



'0 Je 



[ [ W G (s)d£ds) , t e [o,r]. 
Jo ./^ / 

(3.10) 

On the other hand, we recall that condition D(fy) — R is equivalent with 

j* < 00 and lim ^ , = +00. (3-11) 

IpHoo \p\ 

So, there exists N = N(u>) such that 

\ Vx (s)\>N^f( Vx (s))>2C5\ Vx (s)\. 
Consequently, for C > \Qt\ we have 

/ / \vx(s)\d£ds =11 \r) X {s)\dtds+ [ [ MsMds 

Jo Jff J J\v\(s)\>N J J\t) X (s)\<N 

1 rt 




- 2C6 J I f(Vx(s)Hds + NC6. 



e 
12 



Substituting this into (13.101) . since j > 0, we obtain the estimate 

A JO Jiff 

(3.12) 

/ [ (yx + Wc-zxMds^dil + lxW), te[0,T], 
Jo Je 

which implies 

/ /(j(^(s))+j*(r/A(s))^<C 1 (l + |x| 2 _ 1 ), (3.13) 
Jo Je 

and t 

/ f(yx + W G - z x fdids < 2Ad(l + \x\\), (3.14) 
Jo Je 

where C\ is a suitable random constants. 

3.2 Convergence for A — > 

Since (by |HD and 

lim j(m)/|m| = oo, lim j*(w)/|-u| = oo, (3.15) 

\u\— >oo »oo 

we deduce from (13.131) that the sequences {z\} and {r] X } are bounded and 
equi-integrable in L 1 (Qt). Then by the Dunford- Pettis theorem the se- 
quences {z\} and {?7a} are weakly compact in L 1 {Qt). Hence on a sub- 
sequence, again denoted by A, we have 

z\ — > z, Tj\ — ► Tj weakly in L 1 (Qt) as A — > 0. (3.16) 

Moreover, by (I3.12p . (13.141) we see that z = y + Wq where 

y x -> y weakly* in L°°(0,T; if) and weakly in L^Qt). (3.17) 

Note also that by (13.21) we have for every t G [0, T] 

y A (i) - A (J q Ms) + A(j/ A (s) + W G (s)))ds) = x (3.18) 
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and so the sequence 

{J*(Vx(s) + Xyx(s))ds} is bounded in L°°(0,T; H] 
Hence, selecting a further subsequence if necessary (see (13.101) ). we have 



lim f (ri x (s) + \y x (s))ds = [ rj(s)ds weakly* in L°°(0, T; 
Jo Jo 

(3.19) 

So, by fl3~T8|) we find 

y(t)+A I r](s)ds = x a.eJe[0,T]. 
Jo 



(3.20) 



Since 



/ rj(s)ds G C([0, T]; L 1 (^ > )) n L°°(0, T; 
</o 



£ i — ► J^r](s)ds is weakly continuous in H^tff), hence so is t h-> A J^r](s)ds in 
if. So, defining 

y(ty--A [ r](s)ds + x, te[0,T\, (3.21) 
Jo 

y is an if -valued weakly continuous version of y. Furthermore, we claim that 
for A -> 

y A (t) -> y(t) weakly in if, Vie [0, T]. 

Indeed, since ?7a — > rj weakly in L 1 {Qt) and \{y\ + Wg) — > weakly in 
L^iQr) (since it even converges strongly in L 2 {Q T ) to zero by (13.101) ). it 
follows that for every t G [0, T] 

[ (7] X {s) + X{y x (s) + W G {s)))ds -> / 7?(s)ds weakly in L 1 ^). 
Jo Jo 

Hence by (13.181) ) and the definition of fj we obtain that for every t G [0, T] 
(-A)- 1 ^) -> (-A)" 1 ^^) weakly in LH 



Since A > 0, are bounded in H by (13.121) . the above immediately 

implies the claim. 

From now on we always consider this particular version y of y defined in 
(I3.2ip ). For simplicity we denote it again by y\ so we have 

V\{t) -> V(t) weakly in if, V t G [0, T]. 
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We can also rewrite equation (I3.2ip as 

y t (t) - A V (t) = in ®'{Q T ), y(0)=x. (3.22) 
Now we are going to show that 

r](t,0^^(y(t,0 + W G (t,0) a.e. (t,£)eQ T . (3.23) 
For this we shall need the following inequality 

liminf / y x r] X d£dt < / yrjd^dt. (3.24) 

JQt JQt 

To prove this we first recall equation (11.61) which yields 

3x{Vx + W G ) + f x (r] X ) = (y x + W g )t] X , a.e. in Q T , 
and so by (I3.6P and since j% > j*, we have 

j(Vx + W G ) + j*(7] X ) < (y x + W g )t] X a.e. in Q T , 

which yields 

/ (j(yx + W G )+f( Vx )d^dt< [ (y x + W G ) Vx d£dt. 

JQt JQt 

Since the convex functional 



(z,C)-> / (j(z)+f(())dtdt 

JQt 

is lower semicontinuous on L x (Qr) (and consequently weakly lower semicon- 
tinuous on this space) we obtain that 

(j(y + W G ) +j*{r}))d£dt < liminf / y x r} X d£dt + / W GV d£dt. (3.25) 



A->0 



t 



Furthermore, by (13.121) and again by the weak lower semicontinuity of convex 
integrals in L 1 (Qt) it follows that 

j(y + W G ), f(v)£L\Q T ). (3.26) 
15 



On the other hand, since j(u) + j*(p) > up for all u,p G WL (see f ll.7p ). we 
have 

(W G + y)r]<j(y + W G )+f(v) a. e. in Q T . (3.27) 

Moreover, by assumption (12.31) we see that for every M > there exists 
R = R(M) > 0, such that 

j(—y — W G ) < Mj (y + W G ) on Q R 

where 

Q R = {(t,0 e Q T ■■ \y(t,0 + w G (t,£)\ > R}. 

Since j(y + W G ) G L 1 (Qt) we have, by continuity of j, 

]{-y-W G )<h a. e. ing T , (3.28) 

where h G L 1 (Qt)- On the other hand, since j is bounded from below we 
have 

j(-y-W G )eL 1 (Q T ). (3.29) 

Noticing that by virtue of the same inequality (11.71) we have, besides (]3.27|) . 
that 

-(y + W G )r)<j(-y-W G )+f(r)) a. e. in Q T , (3.30) 
by (13.271) and (I3.28P it follows that a. e. in Qt we have 

\{W G + y) V \ < max{j(y + W G ) + f( V ), j(-y - W G ) + f( V )} G L\Q T ) 

and therefore yrj G L 1 (Q T ) as claimed (recall that W G G L°°(Q T )). 
Now we come back to equation (13 .4p which by integration yields 

I (|2/a(T)| 2 _ 1 -|x| 2 1 )+ / y x r) X dtdt + \ [ y x (y x + W G )d£dt = (3.31) 
1 Jq t JQt 

and taking into account that 

y\{T) — > y(T) weakly in H, (3.32) 

we have by (I3.3ip that 

hminf / y xVx dCdt < -\ {\y(T)\\ - \x\\) . (3.33) 
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In order to complete the proof one needs an integration by parts formula 
in equation (13.211) (or (13.221) ) obtained multiplying the equation by y and 
integrating on Qt- Formally this is possible because yr\ £ L 1 (Qt) and 
y(t) £ H~ l (&) for all t £ [0, T]. But, in order to prove it rigorously, one 
must give a sense to (y'(t),y(t)). Lemma [3TT1 below answers this question 
positively and by (I3.33P also proves (13.241) . 

We first note that since j,j* are nonnegative and convex and such that 
j(0) = = j*(0), we have for all measurable / : Qt — > M that for all 
a 6 [0,1], 

3U)£L\Q T )^ 3 (af)eL\Q T ) 

and 

3*{f)eL\Q T )^f{af)eL l {Q T ). 
Furthermore, as in the proof of (13.281) by (12. 3p we obtain 

By ( I2.3P the latter is, however, also true for j*, if / £ L l {Q T ) and a is small 
enough. Indeed by (12.31) there are M, R > such that 

j{-s)<Mj{s) \i\s\>R, 

hence replacing s by (— s) 

^j(s)<j(s) ii\s\>R. 

Now an elementary calculation implies that for all p £ K 

f(-p)<R\p\ + jjf(Mp). 

So 

f(-p/M)<^\p\ + ±j*(p). 
Hence for a : = 1 jM we have 

<f(~af)<^ |/| + ^f(f)eL\Q T ). 

Therefore, y and rj constructed above fulfill all conditions in the following 
lemma since Wq £ L°°(Qt). 
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Lemma 3.1 Let y G C w ([0, T}; H^(ff)) n L l (Q T ) and rj G L\Q T ) 
nL°°{Q,T-H l {ff)) satisfy 

y(t) + A [ r){s)ds = x. (3.34) 
Jo 

Furthermore, assume that for some a > 0, j(ay), j*(arj) G L 1 (Qt)- Then 
yr) G L V {Q T ), 

[ y^dt = -\ (\y{T)W-\xW). (3.35) 

and 

Y £ Z £ -^y V in L\Q T ), 
where Y e , S e are defined in (I3.36P below. 

Proof. We set for e > 

Y £ = (1 + eA)~ m y, S £ = (1 + e^)"^, (3-36) 
where m G N is such that m > max{2, (d + 2)/4}. Then 

Y e G C^T]; n H 2m ~ l {0)) C (T([0,T];i/ V)nC(?)) 

and 

S e G -^ 1 (0, T; W 2,q (&)) 1 Kq< 



d- 1 

Hence F £ S e G L 1 (Q T ) and for e — > 

' Y e {t) -> strongly in H- x {@), V t G [0, T] 

— ► 2/ strongly in L 1 (Qt) 
S e -> 77 strongly in L^Qt) (3.37) 

/ £ £ (s)ds^ / 77(s)ds strongly in Hq{&) V t G [0, T]. 
, jo </o 

We note here that the last fact follows because (I3.34p implies that J' rj(s)ds G 
C w ([0,T\;H%(&)). We have also by (l3~34l) 

Y e (t)+A f ^ e (s)ds = {l + eA)- m x, VtG[0,T], 

JO 
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which implies 

±Y e (t)+AE e (t) = 
at 



and, taking inner product in H l (0) with Y e (t), we obtain 

\Y £ (t)\-i + / Ze(t)Y £ (t)d£ = 0, a.e. t e [0,T]. 



1 d 

2 dt 



Hence 



lim / £ 6 (t)y e (i)d£di = ~\ (Ij/COI-i - \x\\) (3.38) 
£ ^JQt a 

and by (13.371) we may assume that for e — > 

F £ — >■ y, S e -> 77 a. e. in Q T . (3.39) 
Moreover by (ll.7p we have 

a 2 £ £ F e < j(aY £ ) +j*(aZ £ ), -a 2 ^ £ Y £ < j(-aY £ ) + j*(aZ £ ) a. e. in Q T . 

(3.40) 

Now we claim that for e — > 

j(ay £ ) -> j(ay), J*(« s £ ) J*(a»7), j(-«^) -> j{~ a v) in L l {Q T ). 

(3.41) 

By (I3.39P these convergences hold a.e. in Q T . So, in order to prove (13.411) 
it suffices to show that {j(aY £ )}, {j*(aS £ )}, oY £ )} are equi-integrable 
on Qt and so that they are weakly compact in L 1 {Qt)- To this end let 
y G L x {0) and let Y £ := (1 + eA)~ 1 y, i.e. y is the solution to the equation 

Y £ - eAY £ = y, in 0, 

(3.42) 

Y, = 0, on 



It may be represented as 

^(0= / G(£,6)2/(£i)#i, V£g^, (3.43) 



where G is the associated Green function. It is well known that J G(£, £i)d£i 
is the solution to (I3.42p with y = 1 so that by the maximum principle we 
have < J^G(£,6Ri < 1 for all £ e 0. 



19 



We may rewrite Y £ as 

J J 



where 



and so ^(5(^,6)^1 = 1 for all f G ^. 

Then, if j (y) G L 1 (i^ > ) by Jensen's inequality, since j(0) = we have 

j(Y £ (0)< I G(U2)d£ 2 I G(£,6)i(2/(£iMi 



So, we proved that for any y G L l (&) with G L 1 (^ > ) 

Iterating and using the fact that (1 + eA)^ 1 preserves positivity we get for 
all m G N 

+ eA)- m y) < (1 + eA)- m j(y), a.e. in G. (3.44) 

Now let y be as in the assertion of the lemma and Y £ as in (13.361) . Integrating 
over Qt, since (1 + eA)~ m is a contraction on L 1 (^ > ), (I3.44p applied to ay 
implies 

j{aY £ {£,t))d£dt < [ j{ay{&,t))d&dt. 



Taking into account that j(ay) G L 1 (Qt) we infer that {j(aY £ )} is equi- 
integrable on Q T . The same argument applies to {j*(aS £ )}, {j(—aY e )}. 

Then (13.401) implies that sequence {£ £ Y^} is equi-integrable on Qt and 
consequently by the Dunford- Pettis theorem, weakly compact in L 1 (Qt)- 
Since {E e Y^} is a.e. convergent to yrj we infer that for e — > 

E £ Y £ ^yr] strongly in L X {Q T ), (3.45) 

which combined with (I3.38f) implies (13.351) as desired. □ 
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We now prove (13.231) . We have 

J'Oa) - j(u) <rj\(z\-u), V u E R a.e. in Q T . 
Integrating over Qt yields 

/ ](z x )d£dt< / j(u)d£dt + / r] X (z x -u)d^dt, \/ u e L°°(Q T ). 

JQt JQt JQt 

Note that, by the definition of we have 

z x = -\r]x + yx + W G . 

Therefore, since z — y + Wo, by (I3.24p and Fatou's lemma we can let A — > 
to obtain 

f j(z)d£dt- [ j{u)d£dt< [ r]{z - u)d£dt, Vu 6 L°°(Qt)- 

JQt JQt JQt 

Now by Lusin's theorem for each e > there is a compact subset Q e C Qt 
such that (d£ ® dt){Qr \ Qe) < e and y, r\ are continuous on Q e . Let (to,Xo) 
be a Lebesgue point for y, rj and yr\ and let B r be the ball of center (to,^o) 
and radius r. We take 



u, if (t,0 e (Q £ n 5 r ) u (Q r \ Q £ 



Here t> is arbitrary in R. Since it is bounded we can substitute into the above 
inequality to get 



(j(z) - j(v) - r)(z - v))d£dt < / (r](z - v) + j(v) - j{z))d£dt. 
B r nQ e JQt\Q c 

Letting e — > we obtain that 

(j(^) - J'(^) - V{* ~ v))d£dt < 0, V«el,r>0. 
This yields 

j(z(to,x )) <j(v) +r](t ,xo)(z(t ,x ) -v), Vd e R. 
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and therefore r](t ,x ) £ dj(z(to,x )) = ^f(z(t , x )). Since almost all points 
of Qt are Lebesgue points we get (I3.23P as claimed. 

Proof of Theorem 12.31 (Continued). Let us first summarize what we have 
proved for the pair (y,fj) £ L X {Q T ) x L 1 (Q T ). We have 

y £ C w ([0,T];H), [' r,(s)ds £ C w ([0, T]; Hq(0)), 
Jo 

^,()6%(U)) fora.e. £ Q r , 

y{t) + A [ r){s)ds = x, ie[0,T], 
Jo 

j(ay), j*(ay) £ L l {Q T ) for some a £ (0, 1]. 
We claim that (y, rj) is the only such pair. Indeed, if (y, ff) is another then 
j{%{y-y))<\j{ay) + \j{-ay) 

and 

f^(y-y))<lj*(ay) + lf(-ay). 

But as we have explained before Lemma 13.11 the right hand sides are in 
L^^Qt)- Hence y — y, r\ — fj fulfill all conditions of Lemma [37X1 and adopting 
the notation from there we have for e > 

Y £ (t)-Y £ (t) = A /"*(E e (s) — T, £ (s))ds 
Jo 

= f A(E e (s) - E e (s))cfc, te[0,T]. 
Jo 

Differentiating and subsequently taking the inner product in H with Y £ (t) — 
Y £ (t) and integrating again we arrive at 



(1 + eAy m (Y £ (t) - Y £ (t)) = 1 1 [ (Y E (s) - Y e (s))(E £ (s) - t £ (s))d£ds 

o Jg 



1 

2 



1 [ ((l + £ A)- m (y(s)-y(s))(l+eA)- m ( V (s)-f ] (s))dZds J te[0,T]. 
o Je 
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Letting e — > and applying Lemma [3.11 we obtain that for t G [0, T] 

~ \y(t) - y(t)\-i = f I {y{s) - y(s))(v(s) - v(s))d£ds < o 

A JO Jiff 

by the monotonicity of 

Now let us consider the ^-dependence of y and r/. By (13.211) . (13. 23ft we 
know that y = y(t, £, uj) is the solution to equation 



y'(t) - A¥(y(i) + W G (*)(")) = a.e. t G [0,T], 
y(0) = x 

and as seen earlier for t] = t](t, £, a;) as in (13 . 1 6[) 

y eC w ([0,T};H)nL\Q T ), V G L\Q T ) 

i 

ijf^eniO.T];^^)), 



(3.46) 



(3.47) 



and 

V (t,C,u;) eV(y(t,£,u;)) + WQ(t,Z,L;) a.e.(U, W )6Q T xfl. (3.48) 

By the above uniqueness of (y,rj), it follows that for any sequence A — > oo 
we have P-a.s. 

y x (t) - y(t) weakly in H = H~\G), V t G [0, T], 
j/a -> 2/ weakly in L l (Q T ), 

[ r) X (s)ds-* [ r)(s)ds weakly in L x {<ff), Vt G [0,T] 
Jo Jo 

and weakly in Hq(0), a.e. t G [0, T], 

Vx V weakly in L 1 (Qt)- 

Since y and 77 are hence strong L 1 ((5r)-liniits of a sequence of convex conbina- 
tions of y\, t] X respectively, and y\ and r]\ are predictable processes, it follows 
that so are y and r\. In particular, this means that Y(t) = y(t) + Waif) is an 
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if -valued weakly continuous adapted process and that the following equation 
is satisfied 



Y(t)-A [ r]{s)ds = x+ [ G{s)dW{s), te[0,T]. 
Jo Jo 

Equivalently 



dY(t) - AV(Y(t))dt = G(t)dW(t), 



(3.49) 



(3.50) 



Y(0) = x. 

In order to prove that Y is a solution of ( 13.50P in the sense of Definition 
12.11 with G(t) replacing B(X(t)) and to prove uniqueness and some energy 
estimates for solutions to equation (13.501) we need an Ito's formula type result. 
As in the case of Lemma 13.11 the difficulty is that the integral 

V(y)Yd£dt 

might be (in general) not well defined taking into account that ^f(Y),Y £ 
L x (Qt) only. We , however, have 

Lemma 3.2 Let Y as above. Then the following equality holds 
1^)1^ = 1 1*1^ _ f f v (s)Y(s)d^ds 



2 1 v 2 



J6 



(3.51) 



+ /_ {Y(s),G(s)dW(s)}^ + ^ l \\G(s)\\ 2 LHsiL2{nH) ds, P-a.s. 



o 



Furthermore, Y £ CV([0, T]; H) n L 1 ((0, T) x 6 x Q), and rj £ ^(^T) x 
iff x Q) and all conditions (12.61) - (12. 9p are satisfied. 

Proof. By Lemma 13.11 we have that Yrj £ L 1 (Qt)- Next we introduce the 
sequences (see the proof of Lemma 13. 1 f) ) for m £ N 

Y £ = {l + eA)- m Y, S e = (l + eA)- m n. 

For large enough m we can apply Ito's formula to the problem 

dY e (t) + AE E (t) = (1 + eA)~ m GdW{t) 

(3.52) 

y e (0) = (1 +eA)- m x = x e . 
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We have 
1 



Y £ (t)\U= z^l-i 



E £ (s)Y e (s)d£ds 



J0 



+ j\Y £ (s),G £ (s)dW(s)). l + 1 - £ \\G E (s)\\ 2 LHsiLHe))H) ds, te[0,T}. 

(3.53) 

where G £ = (1 + eA)~ m G. Letting e -> (since G L°°(Q T j) we get by 

(E35D 



/ 

■/ Q' 



Yrjd^ds, Pa.s. 



Q 3 



Furthermore 

y e (t) -> y(t) strongly in H~ 1 (iff), Vt G [0,T], 
which by virtue of (13.531) yields (13.51H . if we can show that for t G [0, T] 

P-lim / (Y e (s),G e (s)dW(s)) = / (Y(s),G(s)dW(s)). (3.54) 

We shall even show that this convergence in probability is locally uniform 
in t. We have by a standard consequence of the Burkholder-Davis-Gundy 
inequality for p = 1 (see e.g. [T4"l Corollary D-0.2]) that for Y £ := (1 + 
eA)- 2m Y and 6i,5 2 > 



P 



sup 

te[o,T] 



(y( s ),G( s )dw( s )>- / (y 6 ( s ),G e ( s )(W(s)) 



>5i 



(3.55) 



<^ + P 



l|G( S )f 



L HS {LZ{0),H) 



H) \Y{s)-Y £ {s)Wds>5 2 



Since V G G™([0, T]; if), P-a.s. and (1 + eA) 1 is a contraction on if we have 
sup \Y(s) -F £ (s)|_i < 2 sup |y(s)|^.i, P-a.s.. 

se[0,T] s6[0,T] 

Hence by (I2.14p the second term on the right hand side of (I3.55P converges 
to zero as e — > 0. Taking subsequently 5 2 — > 0, (13.551) implies (13.541) . We 
emphasize that, since the left hand size of (I3.5ip is not continuous P-a.s. 
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(though all terms on the right hand side are), the P-zero set of uo G Q for 
which (13.511) does not hold might depend on t. 

Next we would like to take expectation in (13.511) . Note that because 
|^(t)|?_i is not P-a.s. continuous in t we cannot use stopping times to argue 
that (I3.5ip holds with expectation taken for every summand and the local 
martingale term dropped. We need a more delicate argument here. To this 
end first note that by (I3.48P and (11.61) we have 



rj(s)Y(s)=j(Y(s))+j*(r l (s))>0 1 
hence (13.511) implies that for every t G [0,T] 



(3.56) 



\Y(t)t 1 <\x\\ + N t + t \\G(s)\\ 2 LHs{LHnH) ds, P-a.s., (3.57) 

Jo 



where 



N t := I (Y(s),G(s)dW(s))- 1 , t>0, 
'o 



is a continuous local martingale such that 

< N > t = 2 f \G*(s)Y(s)\l m ds, t > 0, 
Jo 

where G*(s) is the adjoint of G(s) : L 2 (ff) — > H. We shall prove that 



E 



sup |iV t | 
te[o,T] 



< +oo. 



(3.58) 



By the Burkholder-Davis-Gundy inequality for p = 1 applied to the stopping 
times 

t n := inf{t >0:\N t \>N}AT, N G N, 

i/2- 

'o 



E 



sup |iV t | 

t£[0,T N ] 



< 3E 



sup |y(s)UU \\G(s)\\i Hs{L2{ff) . H) ds 

se[o,r N ] 



< 6C E 



sup lYis)^ 
se[o,r N ] 



1/2 



(3.59) 
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where 



C :-- 



E 



\\G(s)\\ 2 LHs ( L 2^. H )ds 



1/2 



< oo. 



Since Y G C w ([0, T]; H), we know that s i— > | is lower semicontinuous. 
Therefore by (13.571) 

|2 



SUp |y(s)|_! 

se[o,Tjv] 



sup |V(s)|li < larlli + sup \N S 

sG[0,Tjv]nQ se[0,rRr] 



||G(s)|| r „„tT2ts?\.u\ds, P-a.s.. 



lLi TS (i 2 (^);fl")"' 5 ' 



So (13~59|) implies that for all N G N 





sup \N t \ 


\ <36C 2 




sup |iV 8 | 


+ C 2 




te[o,T N ] 






S£{0,T N ] 





which entails that 



sup E 

JVeN 



sup \N t \ 
te[o,r N ] 



< oo. 



By monotone convergence this implies (I3.58p . since N t has continuous sample 
paths and | T as N — > oo. Now (13.571) implies that also 



E 



sup ir^i 2 .! 

te[o,T] 



< oo. 



(3.60) 



By ( 133HD . (13361) and fl33TD it follows that 

rjY G L 1 ((0,T) x iff x Q). (3.61) 

Hence by (13.561) 

j(Y),f(rj)eL 1 ((O,T)x0xQ) 

and therefore 

y,7/ G L l ((0,T) x x ft). 

Taking expectation in (13.511) we see that t i— > E[|y(t)| 2 1 ] is continuous. Since 
y G C™([0,T];#), P-a.s., (13T60]) then also implies Y G CV([0, T]; H). This 
in turn together with (I3.49P implies that also (12.71) holds. □ 

Now we come back to the proof of Theorem 12.31 noticing that Lemma 
13.21 also implies the uniqueness of the solution Y and estimate (I2.15p . This 
concludes the proof of Theorem 12.31 □ 
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4 Proof of Theorem 



2.2 



Consider the space 

X = \x e C w ([0,T];H) nL 1 ((0,T) x x ft) : X predictable, 

(4.1) 

sup E[e- 2at |X(t)|^] < M 2 , E / j(X(s)d£ds < M 2 \, 
te[o,T] Jq t > 

where a > 0, Mi > and M 2 > will be specified later. 
The space J?T is endowed with the norm 

/ \ 1/2 

||X|| Q = f sup E[e~ 2at \X(t)W 
\*e[o,T] 

Note that Jfc is closed in the norm || ■ || a . Indeed, if X n — > X in || • || a then 
since 

E / j(X n (s))d^ < M 2 , V rt G N, 

(I3.15P implies that 

X n ^X, in L 1 ((0,T) x ^ x ft) 
and by Fatou's Lemma we get 

E / j{X(s))d£ds < M 2 . 



JQt 

as claimed. Now consider the mapping Y : — »■ defined by 

Y = F(X), (4.2) 

where Y e C w ([0, T]; H) D ^((0, T) x x ft) is the solution in the sense of 
Definition 12.11 of the problem 

dY(t) - AV(Y(t))dt = B(X(t))dW(t) in Q T , 
*(Y(t)) = on S T , (4.3) 
Y(0) = x m0. 

We shall prove that for a, Mi, M 2 suitably chosen, T maps J?T into itself and 
it is a contraction in the norm II • |L. 
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By (USD, (I53B and flU) we have 



^(01-1 + 




(j(Y(s))+j*(r)(s)))dtds 



J0 



(Y(s),B(X(s))dW(s)). 



+\ I \\B(X(s))\\l Hs(L2{ ^ H) ds + ~ IxW, te[0,T]. 



By Hypothesis {H 2 ) we have 



- sup E[e- A *|y(t)|i 1 ]+2e- Art E / (jW+J*(lWM<b 



te[o,T] 




o 



1 . 

< - X 

- 2 1 



1 . 

< - \x 
~ 2 1 



2 

-1 



2 

-1 



L 2 

— sup 

2 te[o,T] 



-2at 



E\X{s)\ 2 _ x ds 



2 te[0,T] 



sup / e- 2a{t - s) ^e- 2as \X{s)Wds < - \x\\ 
e[o,T] Jo 2 



L 2 M 2 
4a 



Hence 



and 



T 2 M 2 

sup E[e- 2at \Y(t)\ 2 _ J < — - + \x 



te[o,T] 



2a 



E / (jW*))+J*(^)))R< 



L 2 M X 2 
2a 



+ \ x \ 



_2aT 



Hence for a > L , Mf > 2|a;|i 1 and M 2 > Mfe the operator T maps 
into itself. By a similar computation involving Hypothesis (-^2) we see that 
for Mi, M2 and a suitably chosen 



||yi-y 2 ||«<-5=||*i-*2|| Q 

1 a 



(4.4) 



where y = TX iy i = 1,2. Hence for a suitable a, T is a contraction and so 
equation X = T(X) has a unique solution in T. This completes the proof. □ 
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